from msq_bistable_feedback_switch import *
import sys

#True:Desktop run; False:Server run;
debug = True

#Simulation running time
time = 16000

#Apply experimentally measured noise in calculations
noise = False

#Input vector dosages [ng]
# 0 [A][pmin][TAL-B:KRAB][BFP]
# 1 [A][pmin][TAL-A:VP16]
# 2 [B][pmin][TAL-A:KRAB][mCITRIN]
# 3 [B][pmin][TAL-B:VP16]
# 4 [CMV][PIR][TAL-B:KRAB]
# 5 [CMV][PIR][TAL-A:VP16]
# 6 [CMV][ETR][TAL-A:KRAB]
# 7 [CMV][ETR][TAL-B:VP16]
# 8 [CMV][PIP:KRAB]
# 9 [CMV][E:KRAB]
copies = [0.0 for i in range(10)]
copies[0] = float(sys.argv[1])
copies[1] = float(sys.argv[2])
copies[2] = copies[0]
copies[3] = copies[1]
copies[4] = 20.0
copies[5] = 20.0
copies[6] = 20.0
copies[7] = 20.0
copies[8] = 100.0
copies[9] = 100.0

#Input signals
#signals[0][...] = pristinamycin input times
#signals[1][...] = pristinamycin input durations
#signals[2][...] = erithromycin input times
#signals[3][...] = erithromycin input durations
signals    = [[],[],[],[]]
signals[0] = [5000]
signals[1] = [2000]
signals[2] = [0]
signals[3] = [2000]


msq = MSQfeedbackSwitch(debug, time, noise, copies, signals)
msq.run()
